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Abstract 


We  derive  a  multipole  scattering  solution  for  a  system  resembling  a  simple 
cell.  In  the  model,  a  spherical  cytoplasm  is  surrounded  by  a  concentric  cell 
membrane.  Contained  within  the  cytoplasm  is  a  nonconcentric  spherical 
nucleus.  Because  of  the  nature  of  the  (multipole  expansion)  solution, 
numerical  results  can  be  acquired  quite  rapidly.  We  show  that  the 
resulting  scatter  is  very  sensitive  to  the  system  geometry  and  optical 
properties.  Such  a  solution  may  also  be  used  to  calculate  the  scatter  from 
fluorescing  molecules  within  the  cell. 
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1.  Introduction 


Light  scattering  is  a  sensitive  tool  that  can  be  used  to  determine  the 
geometry  and  composition  of  particle  systems.  This  valuable  tool  has  been 
used  in  the  laboratory  to  characterize  biological  systems  from  bacteria  to 
red  blood  cells  [1-11].  The  polarization  state  of  the  scattered  light 
expressed  by  the  light-scattering  Mueller  matrix  has  been  shown  to  be 
especially  sensitive  to  small  changes  in  the  system  [5-8]. 

Modeling  the  scattering  from  complex  cellular  systems  is  difficult  owing 
to  fheir  complex  cellular  structure,  and  as  a  result  very  few  programs  have 
been  developed  in  this  field.  One  approach  has  been  to  take  advantage  of 
the  cell's  relatively  small  refractive  index  within  the  host  medium  and  to 
use  the  anomalous  diffraction  approximation  [9,10].  Another  approach 
has  been  to  apply  the  finite-difference  time-domain  technique  to  cells 
[12,13].  Although  these  techniques  have  proven  valuable,  both  are  limited. 
The  anomalous  diffraction  approximation  does  not  accurately  predict  the 
backscatter  from  systems.  It  is  in  the  backscatter  that  polarization 
techniques  have  been  shown  to  be  especially  sensitive  to  small  changes  in 
the  system  geometry  and  chemistry.  Although  the  finite-difference 
time-domain  technique  can  model  virtually  any  system,  it  is  extremely 
computer  intensive. 

In  this  report  we  develop  a  model  that  can  be  used  to  approximate  the 
scatter  from  a  cellular  system.  This  technique  uses  a  multipole  expansion 
of  the  fields  internal  and  external  to  the  cell.  The  model  cellular  system, 
shown  in  figure  1,  is  composed  of  a  spherical  cytoplasm  surrounded  by  a 
spherical  cell  membrane.  Contained  within  the  cytoplasm  is  a  spherical 
nucleus,  which  can  be  placed  at  any  location  within  the  cytoplasm. 
Although  this  model  is  also  severely  limited  (since  the  components  of  the 
cell  are  restricted  to  being  spherical),  it  has  some  advantages:  calculations 
are  very  rapid  and  exact.  This  model  can  be  used  for  examining  the  effects 
of  various  system  parameters,  such  as  the  cell  membrane  composition  or 
thickness,  the  nucleus  size  and  position,  the  inclusion  of  organelles  within 
the  cytoplasm,  or  changes  in  the  cytoplasm  chemistry.  The  derivation  of 
this  model  is  similar  to  one  given  in  previous  works  [14,15]  in  which  we 
examined  the  scatter  from  a  randomly  placed  particle  within  a  spherical 
host.  In  the  present  work,  a  concentric  outer  layer  is  added  to  the  host 
(cytoplasm)  to  represent  a  cell  membrane.  We  show  in  the  results  section 
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that  the  inclusion  of  this  concentric  cell  membrane  surrounding  the 
cytoplasm  can  affect  the  polarization  state  of  the  scattered  light,  especially 
in  the  backscatter  region. 

Finally,  we  note  that  since  this  is  a  complete  solution,  both  the  scattered 
and  internal  fields  are  known  for  the  system.  Using  reciprocity  and  the 
solution  for  the  internal  fields,  one  can  calculate  the  scattered  fields  from  a 
fluorescing  molecule  located  at  any  position  within  the  cell  [16]. 
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2.  Solution 


The  geometry  of  the  scattering  system  is  shown  in  figure  1.  A  spherical 
cytoplasm  of  radius  ai  and  complex  refractive  index  mi  is  centered  on  the 
yii  coordinate  system.  An  outer,  concentric  cell  membrane  having 
outer  radius  03  and  complex  refractive  index  m3  is  also  centered  on  the  xi, 
yii  coordinate  system  such  that  ai  <  03.  A  spherical  nucleus  of  radius 
02  and  complex  refractive  index  m2  is  centered  on  the  X2,  y2,  Z2  coordinate 
system  at  a  position  xi  =0,yi=0,zi=d  such  that  oi  -  02  >  |d|.  In  order 
for  the  scattering  geometry  to  be  completely  general,  the  wavevector  of 
the  incident  radiation  is  oriented  at  angle  a  with  respect  to  the  zi  axis.  The 
wavelength  and  wavevector  of  the  plane  wave  in  the  nonabsorbing, 
nonmagnetic  incident  medium  are  A  and  k,  respectively.  The  complex 
wavevector  in  media  of  refractive  index  m^  is  kj.  To  simplify  the 
equations,  we  take  the  permeability  of  the  spheres  and  the  surrounding 
media  to  be  the  same. 

We  can  find  the  scattering  solution  by  simultaneously  satisfying  the 
boundary  conditions  at  each  interface.  We  consider  the  fields  incident  on 
each  cell  component  separately.  These  fields  are  expanded  in  terms  of  the 
vector  spherical  harmonics,  which  have  the  following  form  in  this 
derivation; 


m: 
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nm.j 
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fCTj 


+  h 
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where  the  index  j  corresponds  to  the  coordinate  system  used  {j  =  1, 2), 
Zn\krj)  are  the  spherical  Bessel  functions  of  the  first,  second,  third,  or 


3 


fourth  kind  (p  =  1, 2, 3, 4),  and 


P-(COS0,)  - 


l{2n  +  l)(n  -  r 
2(n  +  m)! 


-P^{cosej), 


where  P^(cos  0j)  are  the  associated  Legendre  polynomials.  We  assume  a 
time  dependence  of  exp(-iu;f) . 


2.1  Outer  Cell-Membrane  Interface 

We  first  examine  the  fields  that  strike  the  outermost  cell-membrane 
interface  {n  =  03).  We  consider  an  arbitrary  field  incident  on  the  system 
that  can  be  expanded  with  the  spherical  Bessel  functions  of  the  first  kind, 
in(fcri):  ^  ^ 

Einc  =  X]  ^ 

n=0  m——n 

Similarly,  the  scattered  electric  field  may  be  expanded  with  the  spherical 
Bessel  functions  of  the  third  kind,  h^n^{kri): 

E  +  (5) 

n=0  m=—n 

The  fields  inside  the  cell  membrane  may  be  expanded  into  incoming  and 
outgoing  spherical  waves  with  spherical  Bessel  functions  of  the  fourth 

kind,  h^\kiri),  and  third  kind,  hh  \kiri): 

n=0  m——n 

The  application  of  boundary  conditions  at  the  outer  cell-membrane 
interface  for  the  above  three  equations  yields  two  sets  of  equations: 


anmkzi>n{kaz)  +  Cnm^S^n  ^(^03)  —  inmk^nH^S^3.)  +  knm.  k^lrHk3a3), 
anmi’ni^as)  +  (^nm^n^\ka3)  =  inm^n^\haz)  +  ^nm^n'^Hhaz), 
bnm'>Pn{kciz)  +  ~  jnmtn  +  Inm  ^nHhaz), 

hnM'nikaz)  +  dnmkz^^^Kkaz)  =  jnmk^'^^Hhaz)  +  lnmkiT\kzaz), 


(7) 

(8) 
(9) 

(10) 
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where  '4’n{i')  and  in\r)  {q  =  1,2)  are  the  Riccati-Bessel  functions,  defined 
by 

V'n(r)  =  rj„(r)  and  (11) 

and  the  primes  denote  derivatives  with  respect  to  the  argument. 

2.2  Cytoplasm-Cell-Membrane  Interface 

We  now  examine  the  fields  that  strike  the  cytoplasm-cell-membrane 
interface  (ri  =  ai).  The  fields  in  the  region  \d\  <ri<ai  may  be  expanded 
into  incoming  and  outgoing  spherical  waves  with  spherical  Bessel 
functions  of  the  fourth  kind,  h^n\kri),  and  third  kind,  h!^\kn)  : 

oo  n 

^Iph  ^  IZ  J  -I-  /nrnN^rn,!  -|-  gnmMnli,l  +  1-  (12) 

n~0m=—n  ’ 

The  application  of  boundary  conditions  at  the  cytoplasm-cell-membrane 
interface  yields  two  sets  of  equations: 


'>'nrnk\^^\k2,ai)  +  knmkl^^\kzai)  =  &nmkz^^\kiai)  +  gnmkz^^n\kiai),  (13) 

‘^nm^n  \kzai)  +  knrn^n'^\kza\)  =  ^nm^n^\k\ai)  +  gj^yn^'^'^\kiai),  (14) 

jnm^n  {ksO-l)  Inm^n'^  (ksUi)  =  fnm^n^  (kiCli) -j- hnm^n\kiai) ,  (15) 

^nm  ki^n  \kzai)  —  fnm.kz^n^'^ {k\ai) hjimkz^n'^^ {k-yai).  (16) 

Since  our  primary  concern  is  with  the  scattered  fields,  we  can  write  the 
scattered  and  internal  field  coefficients  directly  in  terms  of  the  cytoplasm 
internal  field  coefficients: 

O-nm-^h  ^  ^  'I'  9nmG^\  (17) 

+  dnmDlf^  —  +  hnmHlf\  (18) 

where  J  =  1  or  2.  The  coefficients  can  be  found  from  equations  (7)  to  (10) 
and  (13)  to  (16);  by  applying  the  Wronskian  formula  for  Riccati-Bessel 
functions  [17] 

=-2j,  (19) 
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we  can  derive  the  following  expressions  for  these  coefficients: 

=  A:ife3^„(fca3)^y^(fc3a3)  -  (20) 

=  kk\'il)n{kaz)in^\^^az)  —  kikzi>'ni}^az)^n\kzaz) ,  (21) 

=  kik3^^\kaz)^n^\kza3)  —  kki^J;  \kaz)^li  ^{kzas),  (22) 

=  kki^l^\ka3)^'y\ha3)  -  kikz^'^^^kazj^i^^kzaz),  (23) 

Elf'>  =  kk3^l^\kiai)^'y\k3ai)-kki^'J;^\kiai)(i^'>{kzai),  (24) 

irW  ^  fcfci4^)(/ciai)^;(‘^n^'3ai)-fc^3d^Hfci«i)^i'^^(^3ai),  (25) 

GJ/)  =  fcfe342n^'i«i)d‘^^(^3ai)-fcfcid^Hfci«i)^n'^H^3ai),  (26) 

=  A:fci4^)(fciai)^^^‘^H^3ai)  -  A:A:3d^HA:iai)^^'^^(A:30i)-  (27) 

2.3  Cytoplasm-Nucleus  Interface 

We  now  examine  the  fields  that  strike  the  cytoplasm-nucleus  interface.  We 
examine  these  fields  in  the  0:2,  y2,  ^2  coordinate  system  (j  =  2).  The  fields 
inside  the  nucleus  may  be  expressed  by  the  spherical  Bessel  functions  of 
the  first  kind,  jn(fc2r-2): 

00  n  .  . 

Efnt  =  X]  ^  PnmMnm,2  +  Qnm^nm,2- 

n=0  m~—n 

The  fields  in  the  cytoplasm  may  be  expressed  into  incoming  and  outgoing 
spherical  waves  by  spherical  Bessel  functions  of  the  fourth  kind, 

hn\k\r2),  and  third  kind,  hn\kir2): 

EL=f:  E  rnmMll2  +  wNS,2  +  WMit,2  +  nnmNlti,^  (29) 

n=0  m=—n 

Applying  boundary  conditions  at  the  inclusion  sphere  interface  yields  two 
sets  of  equations: 

Pnmk\'^n{k2(^2^  'f'nm  k2^n\k\a2)  +  tnmk2^nKkia2),  (30) 
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(31) 


Pnm‘^ni^‘iO'2)  —  +  tnm^n^^(fcia2), 

qnmi’n{ha2)  =  Snm^l^Hkia2)  +  Unm^i^\kia2) ,  (32) 

Qnmkl'll>ji{k2a,2)  =  Snmk2^n^\kia2) -{- Uimnk2^n^\kia2) .  (33) 

We  can  eliminate  the  nucleus  field  coefficients  (pnm  and  g„^)  to  find 
relationships  for  the  cytoplasm  field  coefficients.  After  a  little  bit  of 
algebra,  we  have 


kl^n^\kia2)lpnik2a2)  -  k2^n\kia2)lpn(k2a2) 

mn  —  ^nm  7T\ - 77T\ -  = 

k2^n  {kia2)i>n{k2a2)  —  ki^n  {kia2)'ll)n{k2a2) 

_  k2^n^\kia2)lpn{k2a2)  -  ki^n\kia2)lp!n{k2a2) 
k\^n  \kia2)'tpn{k2a2)  -  k2^n  \kia2)tpn{k2a2) 


(34) 

(35) 


The  coefficients  (Q'^  and  Q® )  are  similar  to  the  Mie  scattering  coefficients 
[18]. 


2.4  Fields  in  the  Cytoplasm 


The  fields  interior  to  the  cytoplasm  are  expressed  by  equations  (7)  to  (10), 
while  the  fields  exterior  to  the  nucleus  are  expressed  by  equations  (34)  and 
(35).  These  fields  are  identical.  We  can  equate  these  fields  and  express  the 
coefficients  ejjjTj,  fnm^  dnirf!  and  /inm  irt  terms  of  the  coefficients  rnm.i  Snmi 
tnm,  and  Unm  using  the  translation  addition  theorem.  For  a  translation 
along  the  z-axis  with  no  rotation,  the  vector  spherical  harmonics  are 
related  by 


,,  =  E 

n'=0 


^n' 


n'm.l  ’ 


(36) 


N, 


(«)  _  R(n 

nm,2  2^  -^n' 
n'=0 


("."i)TV/r(9) 


(37) 


where  q  denotes  the  order  of  the  spherical  Bessel  functions  {q  =  3, 4).  This 
relationship  is  valid  in  the  region  where  r  >  |d| .  Explicit  expressions  by 
which  these  coefficients  can  be  calculated  are  provided  in  the  appendix. 
Substituting  equations  (36)  and  (37)  into  equation  (29)  yields 


Cnm  —  ^ 


n'=0 


(38) 
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(39) 


CXD 

Jnm-  —  2-^  ^n'm^n  ^  m^n  ’ 

r?/=0 

oo 

J„  =  “"d  (40) 

n'=0 

oo 

Km  =  Un'm^^n  +  Km^^ 

n'=0 

Substituting  equations  (38)  and  (41)  into  equations  (17)  and  (18)  yields 

,A\P + C„„cw  =  E  [g!,''’  +  «-4'"]  +  K’  +  <«> 


n'=0 

The  exterior  field  coefficients  of  the  nucleus,  tnm  ^ud  Unm/  luay  be 
calculated  if  we  eliminate  the  scattering  coefficients  Cnm  ^i^d  Km  ii^ 
equations  (42)  and  (43): 


E  rp{mm,l)  TT\n 


Jn,m,l) 


bnmPn  =  E  *n'rr,T^ +  Un'mK' 


j{n,m,2) 


where 


A  = 


,(»,™.i)  ^  |c.(2)  [gu)  +  Qr,£(.)]  _  ci,»  [Gi?>  +  o;,£;<«] }  , 

|di«  [bo>  +  -  do>  +  q;.fP>]  }  . 
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uM  =  ^l-'.-)  (flP)  [ifO)  +  qj,4»]  _  b(.)  +  (JJ,J.P)]  J .  (51, 

Since  the  incident  field  coefficients  are  known,  equations  (44)  and  (45) 
represent  two  sets  of  simultaneous  equations  that  can  be  solved  through 
matrix  inversion  techniques  for  the  two  sets  of  field  coefficients.  Although 
the  solution  is  general  for  any  incident  field,  we  consider  specifically  the 
case  of  plane-wave  illumination  whose  wavevector  is  oriented  at  angle  a 
with  respect  to  the  z\  axis,  as  shown  in  figure  1.  Mie  scattering  derivations 
typically  restrict  the  plane  wave  so  that  a  =  0.  Since  we  restrict  the 
nucleus  to  be  centered  on  the  z\  axis,  we  must  remove  the  restriction  on 
the  incident  plane  wave  for  the  derivation  to  be  completely  general.  When 
the  plane  wave  is  polarized  perpendicular  to  the  x-z  plane  (transverse 
electric),  the  coefficients  are  found  to  be  [15] 

^TlTtl  ~  ~  ,  \  \  !  (t),  —  Tn.^fTl.  -C  177  -1-  1  i  { nr\ci 


[V(n-m)(n  +  m  +  l)P^+i(cosa) 


-  •^(n-m-M)(n-l-m)P;;"-i(cosQ)  (51 

2i"+2  ^ 

(.  -  (,™  _  +  1)  [  l(n-m  +  l){n-m  +  2) 

■  n(n  + 1)  [V  ~(2»  +  l)(2n  +  3) 

.  I(n  +  m  +  l)(n  +  m  +  2)  . ,  J 

V  {2n+l){2n  +  3)  ^n+i  i^^osa)  (54 

^  mP^  {cos  a) 

n{n  +  l)  sina  ’ 

When  the  plane  wave  is  polarized  in  the  x-z  plane  (transverse  magnetic), 
the  coefficients  are  found  to  be 


G>nm.  —  d. 


TM_.,TB  . 

'n.m.  ^^n.m  d.IlU. 


"nm  ^nm 


nm  ^nm  * 
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2.5 


Scattering  Amplitudes  and  Efficiencies 

We  consider  the  scattering  amplitudes  in  the  far  field,  where  kr\  >  ka. 
The  scattered  fields  in  this  case  are  in  the  6  and  (p  directions.  In  this  limit, 
the  spherical  Hankel  functions  reduce  to  spherical  waves: 

The  scattering  amplitudes  can  be  expressed  in  the  form  of  the  matrix 

J  -ikn  \  S3  52  /  \  ) 

We  solve  the  scattering  amplitude  matrix  elements  by  expanding  the 
scattered  electric  fields  (eq  (4))  in  terms  of  the  vector  wave  functions  and 
then  expanding  the  vector  wave  functions  (eq  (1))  in  terms  of  the 
polarization  directions.  In  the  far  field,  the  r  component  of  the  electric 
field  becomes  negligible  compared  with  the  0  and  (p  components.  After 
some  algebra,  we  have  the  following: 


51  = 

52  = 

53  = 

54  = 


E  E 

n— 0m=— n 


-‘EE 

n^O  m=—n 
00  n 

-^E  E  H) 


tnupi 


n^irrupi 


77—0  m=—n 
00  jn 

E  E 

n=0m=—n 


irriipi 


sin  01 


(60) 

(61) 

(62) 

(63) 


Using  the  following  relationship  for  normalized,  associated  Legendre 
polynomials,  . 

p-"’(cos0i)  =  (-l)'"Pr(cos^i),  (64) 

we  can  derive  the  following  relationships  between  the  scattering 
coefficients: 

cm  =  i-ircz, 


^TM  _  TM 

^nfh  V  nm  •> 
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(65) 


d^-  = 

^nm  \ 


dlM  =  (-ir^: 


rujTM 


where  m  =  -m. 


The  scattering,  extinction,  and  absorption  efficiencies  of  the  system  are 
defined  as  the  cross  sections  per  projected  area  and  may  be  expressed  as 

(kai)^  ^n(n  +  l)  ^  T  c^m  +  'j  ,  (66) 

^  Ln=l  m=—n  ^  / 

_ o  r  oo  n 

^  V  ‘  /  2_^  [^nm^nm  '  ^nm^nm  '  ^nm  ^nm  '^^nm^nm  )  ?  (67) 

^  ln=l  m=~n  ^ 


—  Qext  Qsca^ 


where  a*„  and  6*^  are  the  complex  conjugates  of  a„m  and  bnm, 
respectively.  The  asymmetry  parameter  5  is  a  measure  of  radiation 
pressure  on  the  system  and  is  a  necessary  parameter  used  in  cell 
levitation.  This  quantity  can  be  expressed  as 


Qscaikd) 


+  n(n  +  2). 


- V  m  Re  4- 

^q'\2  Z_/  X^nm^nm  ^  ^nm  ^nm  J 


l{n  —  m  +  l)(n  +  m  +  1) 
(2n  +  3)(2n  +1) 


X  Re  4-  -I-  ,  jTMJTM*  1 

JXt;  [H^nm‘^n41m  +  “nm“n41m  +  (^nm  <^n41m  +  “ram  “n+lmj  • 

Detailed  derivations  for  the  asymmetry  parameter  and  the  efficiencies  are 
given  elsewhere  [19]. 
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Although  equations  (60)  to  (63)  describe  the  electromagnetic  field 
scattered  by  the  cell,  the  meat  of  the  problem  is  solving  the  simultaneous 
equations  described  by  equations  (44)  and  (45).  For  practical  purposes, 
only  a  finite  number  N  of  vector  spherical  harmonics  are  necessary  to 
describe  the  scatter  from  this  system.  In  our  calculations  we  use  the 
criterion  developed  by  Wiscombe  [20]  and  Bohren  and  Huffman  [18]  that 
N  =  x  +  4a;^/3+  2,  where  x  =  Inaz/X.  Higher  order  terms  are 
insignificant,  and  their  corresponding  coefficients  can  safely  be  set  to  zero. 
We  are  left  with  solving  2N  simultaneous  equations.  Solving  this  resulting 
matrix  is  an  process,  so  for  large  size  parameters  the  computation  time 
slows  down  dramatically.  For  example,  the  CPU  time  required  to  calculate 
the  scatter  (1°  scattering  resolution  in  a  plane)  for  an  x  =  25  cell  is 
approximately  28  s  on  a  SUN  4  workstation  and  approximately  175  s  for 
an  X  =  50  cell. 

We  demonstrate  one  application  of  the  model  by  calculating  the 
light-scattering  Mueller  matrix  elements  for  systems  having  slightly 
different  parameters.  Figure  2  shows  these  matrix  elements  for  cells 
having  different  cell-membrane  refractive  indices  (chemistry).  The 
structure  of  the  total  intensity  (element  5ii)  of  the  scattered  light  is 
essentially  the  same  for  all  three  systems.  The  positions  of  the  maxima  and 
minima  are  virtually  the  same,  but  the  amplitudes  of  the  maxima  tend  to 
increase  (especially  in  the  backscatter  region)  as  the  cell-membrane 
refractive  index  is  increased.  The  polarization  Mueller  matrix  elements 
(matrix  elements  5i2,  S33,  and  534  are  shown)  show  a  marked  difference. 
The  positions  of  the  maxima  and  minima  remain  virtually  unchanged  in 
the  forward-scatter  region,  but  in  the  backscatter  these  minima  and 
maxima  shift,  and  the  amplitudes  are  significantly  different.  These  results 
are  consistent  with  previous  experimental  studies,  which  found  the 
backscatter  region  of  the  polarization  matrix  elements  to  be  extremely 
sensitive  to  small  system  changes  [5-8].  Finally,  it  should  be  noted  that  the 
cell  parameters  of  this  figure  are  not  realistic:  the  sizes  are  too  small  and 
the  refractive  indices  are  too  large.  We  chose  these  particular  parameters 
to  illustrate  the  scattering  effects  seen.  These  same  effects  can  also  be  seen 
in  the  scatter  from  cells  having  more  realistic  parameters,  but  because  of 
the  higher  frequency  oscillations  in  the  scattering  signals  (oscillation 


Figure  2.  Light-scattering 
Mueller  matrix  elements  as 
a  function  of  scattering 
angle  at  three  different 
cell-membrane  refractive 
indices.  For  this  system, 

A  =  0.6328  fjm,  ai  =  2.4 
/im,  a2  =  1.0  /im,  az  —  2.5 
//m,  mi  =  1.05, 1712  ==  1-15, 
and  d  =  1.0  //m. 


frequency  increases  with  particle  size),  it  is  much  more  difficult  to  see  any 
trends  in  the  data. 

Figure  3  shows  the  light-scattering  Mueller  matrix  elements  for  cells 
having  different  nucleus  positions.  These  results  look  similar  to  the 
previous  results.  Although  the  total  intensity  varies  more  than  for  the 
cell-membrane  study,  the  polarization  matrix  elements  still  display  greater 
sensitivity,  and  the  backscattered  light  is  most  sensitive  to  changes  in  the 
system  geometry. 

From  figures  2  and  3,  it  is  apparent  that  the  backscatter  is  the  scattering 
region  that  is  most  sensitive  to  small  changes  in  the  particle  geometry  and 
chemistry.  We  now  examine  just  the  backscatter  from  a  larger,  more 
realistic  cell.  The  system  parameters  are  chosen  to  correspond  to  a  10-^m 
cell  contained  in  an  aqueous  external  medium  (m  =  1.33)  illuminated  in 
the  near  irrfrared  (A  =  850  run),  and  are  similar  to  values  used  in  other 
models  [3,11,13].  Figure  4  shows  the  backscatter  Mueller  matrix  for  three 
different  cytoplasm  chemistries.  Even  though  the  cytoplasm  refractive 
index  change  is  only  1  percent  between  runs,  the  backscatter  signals 
change  significantly.  The  maximum  oscillation  frequency,  which  is 
determined  primarily  by  the  size  of  the  cell,  remains  constant  between 
runs,  but  the  positions  of  the  maxima  and  minima  are  shifted  not  only  in 
the  polarization  matrix  elements,  but  in  the  total  intensity  matrix  element 
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Figure  3.  Light-scattering 
Mueller  matrix  elements  as 
a  function  of  scattering 
angle  at  three  different  cell 
nucleus  locations.  For  this 
system,  A  =  0.6328  //m, 

Qi  =  2.4  /xm,  02  =  1-0  /xm, 
03  =  2.5  ^m,  mi  =  1.05, 

7712  =1.15,  and  m3  =  1.20. 


Figure  4.  Light-scattering 
Mueller  matrix  elements  as 
a  function  of  scattering 
angle  for  three  different 
cytoplasm  chemistries.  For 
this  system,  A  =  0.639  /xm, 
Qi  =  4.99  /xm,  02  =  2.0  /xm, 
03  —  5.0  /xm,  7772  =  1.05, 

7773  =  1.1,  and  d  =  1.0  /xm. 


5n  as  well.  Fifty-percent  differences  in  the  percent  polarization  signals  are 
not  uncommon. 
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4.  Conclusion 


We  have  derived  solutions  for  the  total  field  when  an  incident  plane  wave 
strikes  a  cellular  system.  Because  of  the  nature  of  the  solution,  numerical 
results  can  be  calculated  quite  rapidly  compared  to  other  techniques.  The 
derivation  places  no  restrictions  on  the  refractive  indices  of  the  particle 
constituents  or  the  size  of  the  system;  however,  computation  times 
increase  with  cell  size.  The  model  may  be  used  for  examining  the 
sensitivity  of  scattering  to  geometrical  and  chemical  changes  of  the  cell. 
For  more  complicated  cellular  systems  containing  additional  irregularities 
(such  as  organelles),  it  might  be  possible  to  calculate  scattering  and 
absorption  efficiencies  by  incorporating  effective  medium  approximations 
into  the  model  [21,22].  The  model  may  also  be  used  to  calculate  the  scatter 
from  a  fluorescing  molecule  located  at  some  position  within  the  cell.  In 
our  simulations  of  the  light-scattering  Mueller  matrix,  we  find  that  the 
backscatter  is  extremely  sensitive  to  small  changes  in  the  system 
parameters.  We  demonstrate  that  even  a  very  small  (~1  percent)  change  in 
system  parameters  can  have  a  drastic  effect  on  the  scattered  light.  Since  it 
is  extremely  difficult  to  characterize  cells  to  such  an  extent,  because  of 
geometrical,  orientational,  and  chemical  changes  in  time,  it  is  necessary  in 
many  modeling  simulations  to  calculate  average  scattering  properties, 
which  is  computationally  expensive.  This  report  provides  a  technique  that 
can  greatly  facilitate  such  calculations. 
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Appendix.  Translation  Addition  of  Spherical  Harmonics 


Stein  [23]  and  Cruzan  [24]  derived  translation  addition  theorems  for 
vector  spherical  wave  functions.  For  a  translation  along  the  2;-axis  with  no 
rotation,  the  vector  spherical  harmonics  are  related  by 


)  Y'  A(n.m,q)-^{q)  ,  r>{'n;m,q) ^(q) 

”.2  ^n' 


o  =  V  (A 

^^nm.2  ^n'  ^  ^n' 

n'=0 

where  q  denotes  the  order  of  the  spherical  Bessel  functions  {q  =  3, 4).  This 
relationship  is  valid  in  the  region  where  r  >  |d| .  The  translation 

coefficients  ,4^^?’"*’^’  and  can  be  calculated  from  the  scalar 

translation  coefficients 


^(n,m,g)  _  ^{n,m,q)  hd  j  {n'  -  m  +  l){n'  +  m  +  1)  ^(n,m,<7) 

n'  +  iy  (2n'  +  l)(2n'  +  3)  ^”'+i 

_  /  (^  TlZjin'  -j~  m)  ^(n,m,q) 


i'  V  (2n'  +  l)(2n'  -  1) 


W'-l  > 


Ty{n,m,q)  _  ik\Tnd  „{n,m,q) 

n'  -  n'(n'  +  l) 


The  are  scalar  translation  coefficients.  These  can  be  found  via 

recursion  relations  [19]: 


^2n'  +  lj„/(fcid), 


^n'  ’  ~  —y/2n'  +  ljn>{kid), 
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